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We study charge transport in a granular array with high inter-grain conductances. We show 
that the system exhibits a Berezinskii-Kosterlitz-Thouless crossover from the high-temperature 
conducting state into a low-temperature insulating state. The crossover takes place at a critical 
temperature Tbkt oc E c exp{—g}, where E c is the charging energy of a grain and g 3> 1 is the 
dimensionless inter-grain conductance. A uniformly applied gate voltage drives the insulator into a 
conducting charge liquid state followed by an insulating lattice-pinned Wigner crystal state at larger 
values of the gate voltage. Technically, we establish correspondence between the charge and phase 
representations, employing the instanton gas summation in the framework of the phase model. 

PACS numbers: 73.23.-b, 73.23.Hk, 71.45.Lr, 71.30.+h 
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I. INTRODUCTION 

Granular arrays have recently attracted much atten- 
tion as analytically tractable systems to study the inter- 
play of interactions and scatteringi2i2ii£. The advan- 
tage of granular systems is the possibility to separate 
the scattering -induced quantum interference phenomena 
from electron-electron interaction effects. Quantum co- 
herence is relevant as long as the typical dwell-time in a 
grain Tdweii ~ ^/(<7<^) is longer than the dephasing time 
T v . Here g is the dimensionless inter-grain conductance 
measured in units of e 2 /(2Trh), and 5 is the mean level 
spacing of the grains. For two-dimensional disordered 
interacting systems the dephasing time is known to be& 
Tu> ~ hg/T. The quantum interference is thus suppressed 
if t v < Tdweii, which is the case at T > g 2 5. Hereafter 
we assume this condition to be satisfied and consider the 



temperature range g 2 5 < T < E c (where E c 



7(2C) 



is the charging energy of the grains). This allows us to 
focus on the interaction-induced phenomena, while omit- 
ting the interference (incoherent regime). 

Earlier studies of incoherent two-dimensional arrays 
led to conflicting theoretical result qVi 8 . fn Refs. 013 
the low-tempcrature insulating state was found for suf- 
ficiently small inter-grain conductance, g < 1. Upon 
elevating the temperature, the array was shown to un- 
dergo a Berezinskii-Kosterlitz-Thouless (BKT) transi- 
tion into the conducting state. The transition temper- 
ature Tbkt = Tbkt (5) was predicted to go to zero 
at some critical value of conductance g c » 1.8. At 
larger conductances, g > g c , the metallic state was 
claimed to persist down to zero temperature. In Ref. Q 
the same model as in Refs. 0H was studied in the 
regime of large conductances, 5 > 1, Using a pertur- 
bative renormalization group (RG) analysis, the renor- 
malized inter-grain conductance was shown to behave 
as g => g(T) = g - |ln(i? c £/T) (where 2d is the co- 
ordination number of the lattice). This correction is 
essentially similar (and at T « gS crosses ove»2) to 
the interaction-induced Altshuler-AronoA*i°i conductiv- 
ity corrections known for homogeneous disordered sys- 



tems. At T r~j E c g exp{— dg/2], the conductance is renor- 
malized down to g(T) ~ 0(1). Thus, one may expect 
that the system approaches an insulating state at low 
enough temperature even for large bare inter-grain con- 
ductance. Whether such a "high-g" insulator indeed ex- 
ists, and - if so - its nature and relation (if any) to the 
BKT transition, found 7,8 for "low-g" systems, was not 
clarified. 

In the present work we show how to reconcile these 
findings. In particular, we show that there exists a fi- 
nite Tbkt oc exp{— g} at g 1. To this end one needs 
to go beyond the perturbative analysis^ of the model 
considered in Refs. and include non-perturbative 

- instanton - field configurations. A similar program 
was recently carried out for one-dimensional incoher- 
ent arrays^. The conductivity of one-dimensional arrays 
was found to display activated (insulating) behavior with 
the charge gap ~ E c exp{— <?/4}, which is parametrically 
larger than the temperature where the perturbative cor- 
rections become large. It was also shown that the proper 
low-temperature representation of incoherent arrays is 
that of pinned c/iarge-density wave fluctuations. [This 
should be contrasted with the fl uctua ting phase (or volt- 
age) picture employed in Refs. |4l7l l^|.] The activation 
gap corresponds to the energy needed to create a long 
unit-charge soliton. 

In the present paper the charge representation is de- 
rived and analyzed for the two-dimensional setup. We 
find that the charge excitations are localized unit-charge 
two-dimensional solitons. At g ^S> 1 the solitons in- 
teract logarithmically over a large range of distances. 
This leads to a sharp BKT crossoveniiiiS between a 
low-temperature insulating phase with bound charge- 
anti-charge pairs (and an exponentially small number of 
free charges) and a high-temperature conducting phase, 
where the pairs are unbound. The BKT temperature 
Tbkt (9) remains finite (though exponentially small), 
Tbkt(3) ~ E c g exp{— g}, for an arbitrarily high bare 
conductance g. The zero-temperature quantum phase 
transition at g — g c , found in Refs. thus, does not 

exist. Instead, there is a fast but continuous drop of the 
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transition temperature Tbkt(s) in the vicinity of g ~ 1. 

The issue of whether a classical phase transition or a 
crossover occurs as the temperature is lowered, depends 
sensitively on the details of the model. The true BKT 
transition takes place only if the interactions between 
the charged solitons are logarithmic at arbitrarily large 
distances. This is the case for arrays with inter-grain 
capacitances only (in the absence of the grain's self- 
capacitance, no electric field lines can leave the 2d plane). 
In the presence of the self-capacitance, the interaction 
is logarithmic in a wide, but finite range of distances: 
1 < I < exp{g/2} (hereafter distances are measured in 
units of the array lattice spacing). In the latter case, be- 
low the BKT temperature the array's conductivity is not 
zero (in contrast to the former case), but rather exhibits 
the activation behavior, 



c — g exp 



A 

T 



(1) 



where the activation gap is given by A ~ <?Tbkt 3> 
Ibkt- Upon raising the temperature, above the BKT 
temperature the conductivity sharply increases as 



a = gK exp 



-2b 



BKT 



T-T B 



KT 



(2) 



where K and b are non-universal constants of order 
unitjii^. Finally, above the transition region, the con- 
ductivity crosses over to the perturbative prediction a = 
g — \n(gE c /T). We thus have a generally consistent pic- 
ture based on the BKT physics at any value of bare con- 
ductance, g. 

A gate voltage induces a uniform background charge 
q £ [0,1]. However, for small gate voltages the array 
remains in the particle-hole symmetric state with an in- 
teger number of electrons per dot. The transition into a 
non-uniform state (with a non-integer average number of 
electrons per dot) takes place at a critical dimcnsionless 
charge density q* = A/(2E C ). Its physics is similar to 
the transition from the Meissner to the Abrikosov state 
in type II superconductors upon increasing an external 
magnetic field. In our case the role of the magnetic field 
is played by the gate voltage, q, with q* corresponding 
to the lower critical field H c i. The Abrikosov lattice in 
turn corresponds to the 2d Wigner crystal formed by the 
unit-charge solitons. Such a crystal is easily pinned by 
the lattice. Consequently, at low temperatures, the ar- 
ray is in the insulating phase with a residual activation 
conductivity, associated with the thermal creation of de- 
fects. The Wigner crystal melts at a temperature^ of 
the order of Tbkt, leading to a sharp crossover into the 
conducting phase. 

Methodologically, interacting systems may be modeled 
in two alternative ways: in terms of either phase or charge 
degrees of freedom. The two are canonically conjugated 
and the choice between them is a matter of convenience. 
The phase representation is easier to derive microscopi- 
cally starting from the fermionic tunneling Hamiltonian. 



By this reason it was used in the vast majority of works 
on granular systems and quantum dots^j^L^. We found 
it more convenient, however, to work in the charge repre- 
sentation, which is the natural language to describe the 
insulating phase. In the quantum dot context, the charge 
description was introduced in Refs. Here we em- 

ploy its generalization to 2d granular arrays. We intro- 
duce the model in Sec. ITT1 and show that it exhibits the 
BKT crossover in Sec. IIIII Finite gate voltages and the 
pinned Wigner crystal phase are discussed in Sec. llVI To 
facilitate comparison with the body of work on the phase 
representation, we include the proof of equivalence of the 
two models in Appendix IBl 



II. CHARGE REPRESENTATION 

In this section, we introduce the charge representa- 
tion for incoherent (5 — > 0) interacting arrays. To keep 
the presentation compact, we assume all contacts to be 
single-channel, characterized by a reflection amplitude 
r < 1. The generalization to the multichannel case and, 
in particular, to g ^> 1 is discussed at the end of the sec- 
tion and, in more detail, in Appendix lAl while the proof 
of the equivalence of the resulting charge model to the 
more widely used phase model is outlined in Appendix IBl 



A. Single contact 

Consider a point contact between a quantum dot and a 
metallic reservoir. Such a point contact allows for a small 
number of propagating transverse modes which may be 
thought of as one-dimensional electron liquids (with the 
contact situated at the origin, z = 0). Here we consider 
the case of a single propagating mode, deferring the con- 
sideration of multi-mode contacts to Appendix lAl The 
corresponding one-dimensional electron liquid may be 
bosonized in the conventional and described in 

terms of the bosonic field 9(t, z). Its gradient d z 9(r,z) 
has the meaning of a local electron density. As a re- 
sult, the electron number on the dot may be written as 
N = J °° dz d z 9(r, z) = —9(t,0) and, thus, the Coulomb 
energy takes the form (eN) 2 /(2C) = E c 6 2 (t,0). Finally 
the imaginary-time action of the bosonic field reads 

P ( oc 

S[6(t,z)] = Jdrl J dz [(d T 6) 2 + (d z 9) 2 } 

I- co 

+ E c 9 2 (t, 0) - ^ cos[27t<9(t, 0)]| . (3) 

The last term in this expression describes backscattering 
at the point contact with the reflection amplitude r, while 
D is the electronic bandwidth. 

One may integrate out all degrees of freedom with 
z 7^ 0, retaining the field 9(t) = 9(t, 0) only. The corre- 
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sponding action reduces to 

P 

S W = ^ E + ~ — [ dT cos(27r0(r)), 

71 

where oj to = 27rTm, and we have introduced the Mat- 
subara representation though the transformation 9 m = 
dr 9{r)e~ l ^ mT . The dissipative term, 7r|o; m |^, is gen- 
erated as a result of integrating out the continuum spec- 
trum of the degrees of freedom on the dot. Its appearance 
is a consequence of the assumption that the mean level 
spacing is the smallest energy scale in the model, S — ► 0. 



B. 2d array 



We now generalize the single-contact action, Eq. 10} , 
to the 2d array geometry. To this end we introduce the 
vector index 1 to label the grains. We also introduce 
two fields x ,\(t) and 9 v ,\(t) which describe charge trans- 
port from grain 1 in the positive x and y directions, re- 
spectively. In these notations, the instantaneous electron 
density on the grain 1 is given by the lattice divergence 

v • 0i = e x ,i +ex - e xA + e yA+ey - e y . x ( c f. Fig. nj. with 

the backscattering in the contact, characterized by the 
reflection amplitude r, the action reads 



,5' 



E 
i 



-]T(^ m |0; 2 m + £ c (v-£ 



I,m ) 



Dr x ^ 



dr cos(27T0i j i(r)) 



(5) 



i=x,y 



where D is again the bandwidth. As in Eq. Q , the first 
term in the action JSJl describes the dissipative dynam- 
ics originating from integrating out degrees of freedom 
within the grains, the second term is responsible for the 
charging, and the third one describes backscattering in 
the contacts. 
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FIG. 1: Granular array with a square lattice. The massless 
modes rj, explained in the text, correspond to circular currents 
around a plaquette and, therefore, do not contribute to charge 
transport. 

For a square lattice of linear size M, the array con- 
tains M 2 grains and 2M 2 contacts between them. Con- 
sequently the model is written in terms of 2M 2 bosonic 
degrees of freedom, 9i t \{r). In the limit r — > 0, the masses 
of these modes are provided only by the M 2 charging 
terms E C (V ■ 9\ ym ) 2 . Therefore, if the backscattering is 
neglected, only half of the degrees of freedom of the model 
are massive. To see this explicitly, one may rewrite the 
2d vector field 9 through two scalar fields 9 = V^+Vx r\ 
and notice that the charging term contains only the field 
X while the curl-field r\ fully decouples from it. 

In order to find an effective low-energy theory, we shall 
proceed with the renormalization group (RG) scheme 



based on the integration of the high frequency Matsubara 



modes, D' < 



< D, accompanied by the appropriate 



change of the backscattering amplitude r. As long as the 
coefficient in front of the cosine-term in Eq. JSJ) is less 
than the running bandwidth D' , one may treat the fields 
as Gaussian, governed by the first two terms in Eq. J5J. 
As a result, the backscattering amplitude renormalizes 
as 



Dr 



Drexp J v ; 102 



(6) 



where the averaging in (9 2 ) is performed over high- 
frequency fluctuations. Passing to the momentum rep- 
resentation and taking into account both x an( i V com- 
ponents of the fluctuations, one finds 



T 

AM 2 



M 



E E 

\uj m \=D' qn,g y =l 



1 



1 



1 



Eq + ir\LO r , 

1 



7r|w m | 
D 



d n . 

, , In h m — = — 77 In ■ , . 

4tt 2 V E n D' J 2tt 2 JE~U ' 



(7) 



where E^ = AE C J2i=x y sm ( 7r( 7i/(2-^0) is the mass spec- 
trum of the x modes. In the second line we have assumed 
that D' < E c (in the opposite case, D' > E c , one should 
substitute \fE c D' by D'). Notice that the presence of 
the lower limit, D' , in this expression is due to the mass- 
less rotational modes of the field r\. [Note also the dif- 
ference with the ld-system, where all modes are mas- 
sive and, therefore, the result corresponding to Eq. Q 
is independent on the lower limil£.] Combining Eqs. © 
and J7|), one finds that upon integrating out the high- 
frequency modes, the coefficient of the cosine potential 
renormalizes as Dr =>• \JE C D' r. As was discussed above, 
this procedure works as long as \JE C D' r < D' , that is 
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D' > To, where Tq — E c r 2 is the "freezing" temperature. 
For smaller bandwidths, the cosine-term itself provides 
a mass for the rotational modes ?y. As a result, all modes 
acquire a mass and, thus, the renormalizcd backscatter- 
ing amplitude looses its sensitivity to the lower limit D'. 
Therefore, we arrive at the conclusion that for D' < To 
the cosine amplitude saturates at a value about To. In- 
tegrating in this way all Matsubara components, except 
the static one, m = (it is obvious from Eq. J7J that the 
m = component can not be handled in the same way), 
one obtains an effective classical model with the action 

Scl[#1 = |E {( V • ft)' - ^ E COS ( 2 ^.0}' ( 8 ) 
1 i—x 7 y 



where j(T) = 2tt v / T/E c r for T > T and 7 (T) ~ 27rr 2 
for T < T . 

So far we have formulated the model for an array with 
single-channel contacts. Generalization to the N > 2 
channels is achieved by introducing bosonic modes 9i a {r) 
for every channel a = 1 . . . N. One may then integrate 
out the antisymmetric modes (including the spin modes) 
for every contact, retaining only the symmetric (charge) 
mode 6\ — $Z 0=1 see Appendix lAl for details. The 
main result of such a procedure may be summarized 
by the redefinition of the effective backscattering ampli- 
tud o 1 Ti 1 ? r =>■ cm IIa=i r a in the action of the charge 
mode (cat is a numerical coefficient with a finite limit 
Cqo). Consequently the characteristic freezing temper- 
ature changes to To ~ E c Y\. a =i r a ■ Then, the charge 
mode may be described by the same effective classical 
model Eq. © with 7 (T) ~ y/rjE~ c flLi r " for T > T 
and 7 (T)~ 1^=1^ for T<T . _ 

A model which adequately describes an array of metal- 
lic grains assumes that the contacts between grains con- 
sist of a large number N of weakly transmitting chan- 
nels^iiifi. For sufficiently large N, the total conduc- 
tances of the junctions may still be high, g = X^=i > 
1, where t\ — 1 — r 2 <C 1 is the transmission probabil- 
ity in channel a. In this case, one finds that Y[ a r a = 
ex P {I Eq ln (! ~ t2 a )} ~ exp{-.g/2}. Employing the ex- 
pressions derived above, one obtains that the freezing 
temperature is of the order To ~ gE c e~ 9 , while the effec- 
tive amplitude of the cosine potential in Eq. (jSJ is given 
by 

7( t } J^ e ~ 9/2 vl ; t> To , 

[ge-s; T < T . 

In Appendix [B] we show how Eqs. JSJ and (J5J) follow 
from the phase modelSiMi£, demonstrating that the two 
models based on the charge and phase representations, 
respectively, are reduced to the same effective classical 
system. Notice that, since the charge model was derived 
for r„Cl, while t a <C 1 is assumed in the phase model, 
the coincidence of 7 (T) may be expected at best with ex- 
ponential accuracy. The algebraic pre-exponential func- 



tion of g is a result of the evaluation in the framework of 
the phase model, see Appendix IB1 

III. BKT TRANSITION 

In this Section we analyze the physics of the classical 
charge model specified by Eqs. JSJ and Two issues 
are discussed: (i) the spectrum of its charged excitations 
and their interactions, and (ii) the low-frequency charge 
dynamics and the dc conductivity of the array. We finally 
put our findings in perspective by comparing them with 
the results of previous studies. 

A. Charge spectrum 

The lowest energy configuration of the action JSJ) is 
given by 9 = (mod 1) everywhere. Localized exci- 
tations must have integer 9 far away from the core to 
minimize the cosine potential. The total charge of such 
localized excitation is e J (cPl) V ■ 9 = ej ds- 9, where the 
line integral on the r.h.s. is calculated over a distant con- 
tour enclosing the excitation. It is clear therefore that the 
charge of the excitation is quantized in integer numbers 
of e. The simplest (and only stable) charged excitations 
have charge ±e. They consist of a large (i.e., spread 
out over ~ l/ 7 3> 1 grains) localized 2d soliton of unit 
charge, connected to a Id string of links with 9 L = 1. 
The other end of the string may either go to the system 
boundary or be terminated by an anti-soliton with charge 
— e. The soliton solution centered at 1 = can be written 
in the form 9\ = 1 — 1?(1) for the links along the string 
and 9\ = $(1) everywhere else, where |#(|1| — > oo)| — > 0. 
Minimizing the action, Eq. (JSJ, with respect to §, one 
finds the saddle point equation for the soliton solution, 

V (V • $) - ^- V Bin(27n?i)ei = 0. (10) 

i—x.y 

Except for a domain consisting of 0(1) links closest to 
the core of the soliton, ■& is small, justifying an expansion 
of the sine-term in the saddle point equation. As a result, 
Eq. (fTU|l takes the form V(V-i?) — 7 # = 0. Its unit-charge 
solution is: 

*» - 

where K\ is a modified Besscl function, and £ s = 1/^/7 
1, justifying the continuum approximation; finally e; = 
1/Z. The solution is normalized as J (d 2 l) V-0 = 1 to obey 
the charge quantization. 

Substituting this solution back into the action, Eq. JSJ), 
one finds that the soliton energy originates primarily 
from the cosine potential part of the action and is given 
by A ~ E c [y(T)/2n] ln£ s . The large logarithmic factor 
ln£ s = — i In 7 ^> 1 is due to the oc 1/Z 2 behavior of the 
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charge density in the wide range of distances 1 < I < £ s . 
At larger distances, I > £ s , the charge density decays 
exponentially. As a result, the solitons interact logarith- 
mically up to a distance £ s beyond which the interaction 
is exponentially screened. Since the density of thermally- 
excited solitons is n s sw exp{— A/T}, the mean distance 

— 1/2 

between them is l s — n s rts exp{A/(2T)}. It becomes 
comparable to £ s at T w A/(21n£ s ) = E c j(T)/(4w). 
This condition is satisfied at temperatures about the 
"freezing" temperature, T ~ To. Thus, at T < To, the 
thermally-excited charges are essentially non-interacting, 
while, at T > To, there is a neutral (in average) gas of 
logarithmically interacting solitons and anti-solitons. 

In the latter regime the partition function of the 
charged degrees of freedom may be written therefore as 



phase. 



* = £{r/(^). 



n=0 



(d 2 l r , 



■ 2nT 



k.k' 



-VI 



(12) 

where In/ ~ E c -f(T)/T is the fugacity of the logarith- 
mic gas, originating from the solitons core energy. The 
plus/minus signs in the exponent correspond to soliton 
soliton and soliton-anti-soliton interactions, respectively. 

It is well known that the Coulomb gas in 2d described 
by Eq. (|T2)l undergoes the BKT transition 11 ! 12 at a crit- 
ical temperature Tbkt ~ E c ^ / (Air). For T < Tbkt, 
the charges are bound in charge-anti-charge pairs. The 
residual density of free charges is exponentially small and 
given by n s ss exp{— A/T}, where A = Tbkt^^ s 2 . The 
value of A is finite but large, as long as the solitons 
interact with each other logarithmically over a broad 
range of distances £ s >■ 1. Notice that the Coulomb 
interactions in our model are strictly on-site (only the 
self-capacitance, C, is included). The long range of the 
soliton-soliton interactions is due to the fact that in a 
strongly coupled array, j > 1, the charge is spread over 
a large distance £ s ~ exp{g/2}. 

One can modify the model to include mutual capac- 
itances C between neighboring grains (and thus to in- 
clude long-range Coulomb interactions). It is straight- 
forward to show that such modification alters the range 
of logarithmic interactions as £ s — > £ B yl + C'/C, while 
the charging energy now reads E c = e 2 /(2(C + C')). In 
the limit C — ► 0, while C remains finite, the interaction 
range diverges, £ s — * oo. In fact, this was to be an- 
ticipated: since without the self-capacitance no electric 
field lines can leave the system, one deals with the true 2d 
Coulomb interaction, which is logarithmic. In this case 
A — > oo and the density of free charges below Tbkt is 
strictly zero. This is the case of the genuine BKT phase 
transition. For non-vanishing self-capacitance, C > 0, 
the interactions are screened at distances exceeding £ s . 
Therefore, the density of free charges is finite at any 
temperature and the phase transition is smeared into a 
sharp crossover. Above the transition/crossover temper- 
ature the densit y of free charges r apidly increases as^i 
n s ~ exp{— 2by^T-Q^/(T — Tbkt)}, where b is a con- 
stant of order unity, driving the array into the conducting 



B. dc conductivity 



In order to discuss the dc conductivity of the array, 
one needs to restore the low frequency, u -C T, dynamics 
of the classical charge model, Eq. © . This may be done 
formally by keeping the dissipative dynamical term in 
the action. Notice that in the multichannel case, see Ap- 
pendix 1X1 the coefficient in front of \u> m \ acquires a factor 
A" 1 , where N is the number of channels. In the presence 
of strong backscattering, it actually reads Trg~ 1 \uj m \9 l 2 m 
and corresponds to the conventional Ohmic dissipation. 
Since we focus on the low frequencies, it is convenient 
to pass to the Keldysh representation (to avoid dealing 
with the analytical continuation) and consider its semi- 
classical limit. The latter is known to be equivalent to a 
certain Langevin dynamics 19 . 

Here we prefer to take a more phenomenological 
route, leading to the same conclusions. Let us consider 
the static equations of motion following from Eq. JSJ: 



2ttC 



sin(27r6» 4! i) = 0. Since • 0i = V\ is 



the voltage on grain 1, the equation simply expresses the 
fact that in the absence of charge quantization, 7 — > 0, all 
grains are equipotential: VVi = V\ +e .—V\ = 0. Once cur- 
rents are allowed to flow in the array this condition should 
be substituted by the Kirchhoff law, V\+ e . — V\ = RI^i, 
where R = 2nh/ (e 2 g) is the contact resistance, and 
T,i = edt&i,i is the current flowing between grains 1 and 
1 + e,-. Restoring also the 7-term in the equation of mo- 
tion, one thus finds 



• d.0 - E c 



V(V-#)- ^ VeiSin(27r^) 

Z7T • 



■E + £(i). 



(13) 



On the right hand side we have included an external elec- 
tric field E, as well as the Gaussian noise, £(i), with the 
correlator 

27rT 

(*)&',!'(*')) = S(t- t'^^i, , (14) 

9 

to satisfy the fluctuation-dissipation theorem. 

Our goal is to calculate the current, /, in presence 
of a weak uniform field, E. To this end we employ 
Drude-type arguments, saying that I = en s v, where n s 
is the carrier concentration and v is their drift velocity. 
The only mobile carriers in the system are the solitons, 
Eq. (|11|) . whose concentration, n s , we have discussed in 
detail above. Now we concentrate on the drift velocity, 
v. We look for a solution of Eq. I|13f) (without the noise) 
in the form 0(1, t) = 6 (l-vt) + 0*i(l-vi) + a. Here 
6*o (1) is the static soliton solution in the absence of the 
external field, whereas 61 ~ E is a small modification 
of the soliton's shape due to the presence of the exter- 
nal field. Finally the constant vector a is determined 
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by the shift of the minimum of the periodic potential in 
the field: E c -f sin(27rai) = 7rE. Choosing E = Ee x 

and v = ve x , and linearizing Eq. (|13() . one finds that Q\ 
satisfies the equation 

— v — E 
E cF<ff nl 0i = -d x 9 sin 2 (7r6'o,x)e s , (15) 

\"0i g n 

where F { § a }0i = V(V-0i) - 7 ^ e i cos(27r6> , 4 )6> M . The 
velocity, u, is determined by the condition that the r.h.s. 
of Eq. (|15|1 is orthogonal to the translational zero-mode 
of the operator fig \i given by d x 0Q. This requirement 
leads to v ~ g-E. Finally, the dc conductivity is given by 
cr ~ gn s (T). 

As a result, all the conclusions, drawn above, regarding 
the BKT transition/crossover in the soliton density may 
be directly translated to the array's conductivity. In par- 
ticular, for the self-capacitance model at T < Tbkt (em- 
ploying that at low temperature In 7 ~ g) we find Eq. IjTjI. 
i.e. a ~ gexp{— A/T}. Above Tbkt, the conductiv- 
ity behaves as a ~ gexpj— 2&-\/Tbkt/(T — Tbkt)}, see 
Eq. (O. At even higher temperatures, this behavior 
crosses over to the result of the perturbative calcula- 
tion, a = g — \n{gE c /T). 

C. Phase diagram 

Using the results of the previous sections, we are now in 
a position to discuss the phase diagram of the array. An 
array having inter-grain capacitances C only exhibits a 
BKT phase transition between the low-temperature insu- 
lating and the high-temperature conducting phases. Its 
phase diagram on the plane temperature vs. bare inter- 
grain conductance, g, is shown in Fig. [21 Unlike previ- 
ous works« that predicted a zero-temperature metal for 
g > g c ~ 1, we find that the low-temperature phase is 
an insulator for arbitrarily large g. The critical temper- 
ature, Tbkt (ff), however, drops sharply at g ~ 1 and, at 
large g 3> 1, behaves as Tbkt ~ E c gexp{—g}. As shown 
in Appendix [5] the disagreement is not a consequence 
of the different model we use, but can be traced back to 
the disregard of the quantum fluctuations of phase in the 
earlier works. By contrast, Ref. uses a perturbative 
renormalization scheme that neglects instanton configu- 
rations. However, it is precisely these instanton config- 
urations that reflect the discreetness of charge which is 
the key point in identifying the transition. 

In the presence of t he self-cap acitance, C > 0, the 
screening length £ s ~ \/l + C'/C exp{g/2} is finite, and 
the transition is smeared into a crossover. The crossover 
is sharp as long as £ s ^> 1. Regardless of the ratio C/C, 
this is the case for g ^> 1. In this regime the charge 
gap is parametrically larger than the crossover temper- 
ature, A ~ J/Tekt, and therefore the residual conduc- 
tivity below the crossover, though finite, is exponentially 
small, a < gexp{— g}. As a result, quantitatively, there 
is little difference between the models with and without 




FIG. 2: BKT temperature as a function of g. For j > 1, 
Tbkt is exponentially small, but remains finite, a) C = 0: 
the true transition exists for any value of g. Conductance is 
zero below the transition temperature, b) C > 0: crossover 
takes places. It is sharp only if g ^> 1 and/or C' 3> C. In the 
regime T < Tbkt, the system shows activation behavior with 
the gap A = Tbkt m(£ s 2 ). 



self-capacitance. This is not the case for g < 1: unless 
C C, the crossover is gone, and the conductivity fol- 
lows the simple activation law a ~ <?exp{— A/T} with 
A ~ E c . 



IV. FINITE GATE VOLTAGE 

So far we have restricted ourselves to the case of zero 
gate voltage only. A finite gate voltage induces a contin- 
uous background charge q cx Ugatc on the grains. In this 
case the charging term in the action Eq. JSJ has to be 

replaced with S^&q] = E c /TJ2i(V ■ #i - q? ■ Alter- 
natively one may shift the 9 field by q\ to move the q— 
dependence into the pinning term cos(27r(6>i + qh)). 
Since grain coordinates U take only integer values, the 
model is periodic in the g-space with unit period. In 
this work we restrict ourselves to a uniform gate volt- 
age q(\) = q = const, leaving considerations of a random 
background charge for future studies. 

For small q, the system is in a particle-hole symmetric 
"neutral" state: the ground state is still (as for q = 0) 
characterized by 9\ — 0. At some finite value q = q* , a 
transition takes place, where the ground state becomes 
charged (with a non-integer average number of electrons 
per dot) and spatially non-uniform. To find q* , let us 
compute the soliton energy in the presence of q. Since 
the g-dependence of the Hamiltonian is a pure boundary 
effect, one immediately finds the soliton energy A(g) = 
A(0) - 2qE c . At q* = A(Q)/(2E C ) the soliton energy 
A(q) vanishes. This marks the transition into the charged 
state: for q > q* , solitons are created at no cost. 

In the charged phase, the density of solitons in the 
array is finite even at zero temperature. In order to find 
the soliton density at q > q* , one has to take into account 
interaction between the solitons. At small densities, n s < 
= 7, the interaction between solitons is exponentially 
weak. The energy cost associated with a soliton density 



n s reads 



array 



E < (n s ) = n s A(q) + —E c "fn s K 



(16) 



where the second contribution is given by the interac- 
tion energy of a pair of solitons separated by the dis- 
tance l/y/n s , multiplied by the soliton density. Since 
yjj/n s 3> 1, we can use the asymptotic expression for 
the Bessel function, Kq{x) ~ a; -1 / 2 exp{— x} for x — ► oo. 
The optimal density is determined by the minimum of 
E < (n s ). Minimization of B<(n s ) with respect to n s 
yields n s (q) ~ 7/ In 2 [7/(9 — 9*)L where we used that 
A(q) = 2E c (q* — q). Thus, at q > 5*, the soliton density 
rises rapidly until at q — q* ~ 7, the distance between 
solitons reaches £ s . For n s > £,r 2 , the solitons start to in- 
teract logarithmically. Consequently, the expression for 
the energy has to be modified as 



00 

£>(n 8 ) = n s A(q) + ^-E cl n 2 s J Idl K (^yl) 

~ n s A(q) + -^ E cn a (n s - 7 ln^^^ , (17) 



where the second contribution describes the interaction 
energy of the solitons with density n s ; in the volume £ 2 = 
I/7, the interaction is logarithmic [Kq(x) ~ — lnx for 
1 « 1]. In this regime, the minimization yields n s (q) ~ 
2iT(q -q*) + (7/4) ln[(q - q*)/ 7 ] for q - q* » 7. 

Naively, one would expect the system to be no longer 
insulating once the density of solitons becomes finite at 
q > q* - which would be the case if the solitons were 
mobile. However, even though the soliton density in the 
system is finite, n s > 0, it turns out that - except for 
a narrow region q — q* < 7, where the interaction be- 
tween solitons is exponentially weak - the solitons form 
a Wigner crystal which is pinned due to the underlying 
lattice structure. Thus, transport is still activated. 

To understand this fact, we use the analogy with 
the formation of vortices in a type II superconducting 
film 2 ^. The field 9 may be viewed as A x n z , where A 
is the vector potential and n z is a unit vector normal 
to the film. Since the local magnetic field is given as 
h = /in 2 =Vx A, the correspondence goes asV-9 = h 
and the charge quantization in the array is equiva- 
lent of the flux quantization in the superconductor, 
J (d 2 l) h — k (k 6 Z), where h is measured in units of 
the flux quantum <pQ. In this analogy, the gate voltage 
translates to the external magnetic field, H, and the 
gate voltage q* corresponds to the critical magnetic field 
H c i, where it becomes energetically favorable to create 
vortices. The correspondences are summarized in the 
following "dictionary" : 



charge V ■ 6 

& = VV7 

background charge q 

q 



superconducting film 



A x n z = -A 2 V/i 
local magnetic field h 
penetration depth A 
external magnetic field H 

Above H c i there is a finite density of vortices in the 
system which at low enough temperatures form an 
Abrikosov lattice 21 . In a clean film, the vortex lattice is 
free to move, but it is easily pinned by the system bound- 
aries, the underlying lattice structure (as in Josephson 
junction arrays) or any sort of disorder. Upon increasing 
the temperature the vortex lattice eventually meltsii*i^, 
and above the melting temperature T m most of the vor- 
tices are free to move. The melting temperature at finite 
H > H c i is smaller than, but parametrically the same as 
the Berezinskii-Kosterlitz-Thouless temperature at zero 
magnetic field 14 . Thus, at T < T m the system is super- 
conducting while at T > T m the moving vortices lead to 
dissipation. 



charge 
solid 




T m T BKT 



FIG. 3: System properties as a function of an external gate 
voltage. The system becomes charged at q > q* . If the soliton 
density is larger than l/£ s 2 , charges arrange into a (pinned) 
Wigner crystal, a) Phase diagram, b) Conductivity at T < 
T m as a function of q. 

Translating back to our problem this means that at 
q > q* the solitons form a Wigner crystal once their den- 
sity is sufficiently large such that the interaction is loga- 
rithmic. Only in the narrow interval q* < q < q* + £,T 2 
the system is in the conducting charge liquid state. Upon 
increasing the gate voltage, the Wigner crystal forms and, 
due to lattice pinning, the system is an insulator at tem- 
peratures smaller than the melting temperature. The 
latter is of the order of Tbkt- Note that while for q < q* 
charge is carried by individual (thermally-activated) soli- 
tons, for q > q* the mobile charges are lattice defects, 
whose core energy is proportional to the logarithm of the 
lattice constant of the Wigner crystal. 



V. CONCLUSIONS 

The conductivity of a granular material with small 
(g -C 1) inter-grain conductances is controlled by the 
Coulomb blockade effect in separate grains^. Charge 
transport in such an array occurs by electron hops be- 
tween single grains. Because of the Coulomb blockade, 
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the granular array behaves as an insulator at low tem- 
peratures. The characteristic energy determining the ac- 
tivation of the charge transport, is associated with the 
single-grain charging energy E c . 

In this paper, we investigated the properties of a gran- 
ular array at high inter-grain conductance, g 1. We 
concentrated on the simplest model, neglecting the spac- 
ing between the discrete electron levels in the grains 
(8 — ► 0), and in the main part of the paper we also as- 
sumed the ideal limit of zero background charge (<j> = 0). 
We found that the granular array at sufficiently low tem- 
peratures remains an insulator even in the case of g 3> 1 . 
The large inter-grain conductance, however, does affect 
the nature of the charge carrier. Instead of an integer (in 
the units of e) uncompensated charge sitting on a sin- 
gle grain, it is rather a charge-e soliton involving many 
grains. There is a sharp crossover to low conductance 
at temperature the Tbkt oc i? c cxp{— g}. Below the 
crossover, the electron transport is associated with the 
activation of solitons with charges ±e; the corresponding 
activation energy is relatively high, A ~ ^Tbkt- The ap- 
proach to the crossover region from the high temperature 
side can be described in terms of the correlation ra- 
dius for the Berezinskii-Kosterlitz-Thouless transition, 
see Eq. @. Comparison of Eqs. and shows that 
the crossover width is ST ~ Tbkt/s 2 <C Tbkt- 

The effect of a gate voltage q applied uniformly 
throughout the array is to much extent similar to that 
of an external magnetic field applied to a type II su- 
perconductor. Until q reaches a certain critical value 
q* the system's behavior does not change qualitatively 
(that is, it exhibits the BKT crossover from insulator to 
metal). The charge gap, A, and the crossover tempera- 
ture, Tbkt, decrease with increasing q. At q — q* the 
gap vanishes (while the crossover temperature remains 
finite), and for even larger gate voltages there is a cer- 
tain ground-state density, n s = n s (q), of charge solitons. 
As long as this density is small n s (q) < £~ the solitons 
are in a liquid state and the array conducts. At larger 
density, n s (q) > £r 2 , the solitons form a Wigner crystal 
pinned by the lattice. As a result, the array is again in 
the insulating state with the charge gap determined by 
the cost of a defect in the Wigner crystal. 

Finally, let us mention related issues that are not ad- 
dressed in the present paper. The first one is the role of 
disorder. The most relevant is charge disorder equivalent 
to a grain-dependent gate voltage q\. In the extreme sce- 
nario one may assume that q\ are independent random 
variables, uniformly distributed in [0, 1]. One would then 
like to solve the classical statistical problem formulated 
by Eqs. (|13l) and (|14|) with random q\ in the argument 
of sin(27r(6'i ! i + q\)). Despite many similarities with the 
vortex physics, one can not simply transfer the known 
results from the pinned vortex lattice literature^ 2 . The 
reason is that random charge, q\, translates into a strong 
(of order oiH^) fluctuating magnetic field, rather a than 
fluctuating pinning potential. 

Another unaddressed issue is the role of quantum co- 



herence, which enters the problem through the mean level 
spacing S. Our results are valid as long as g 2 S < Tbkt- In 
the opposite limit, the quantum coherence effects (most 
notably Anderson localization) start to interfere with the 
effects of electron-electron interactions, considered here. 
One may expect that both effects drive the system to- 
wards the insulating ground-state. (It is worth mention- 
ing that in both cases the characteristic length happens 
to be exponentially large in the bare conductance, g.) 
The structure of such insulator is not known currently. 
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APPENDIX A: MULTI-CHANNEL CONTACTS 

In the main text (Sec. Ill B|) . we derived a classical 
model for the single-channel case. Here we discuss its 
generalization to N > 2 channels. For every channel 
a = 1, . . . , N of a multi-channel contact, one introduces 
a field 9 a (r). Consider an M x M array with N channels 
in each of the 2M 2 contacts. The quadratic part of the 
action reads 

52 = ^ E M^i E +^(E V - *U) j >( A1 ) 

l.m \ a a / 

while the backscattering is described by 

P 

E f 

S r = ^ E rila / dT COS ( 27r #U,a), ( A2 ) 

I, a i=x,y o 

where the high-energy modes E c < \u m \ < D have al- 
ready been integrated out. [At energies larger than E c , 
all modes are decoupled and, thus, can be integrated out 
for each channel separately.] Here, in order to clarify the 
following evaluation scheme, the reflection coefficients r 
have been given indices specifying the direction, contact, 
and channel. 

Only the 2M 2 symmetric modes 9\ = J2 a #i jQ cou- 
ple to external parameters, such as gate voltages. We 
want, thus, to find an effective action for 9\ by integrat- 
ing out 2M 2 (N — 1) asymmetric modes. To this end 
let us change variables from 9\^ a (a = 1 . . . N) to 9\ and 

9l,a = *l,a-(01-£a'>« WO^ 1 ) («,«' = !■■• N-l). 

While the symmetric fields 9\ are massive due to the 
charging term, all the asymmetric fields 9\. a are massless. 
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As a result, the perturbation theory in powers of rn a con- 
tains only the terms that do not have massless fields B\ a 
in the exponents (cosines). Rewriting the backscatter- 
ing action in terms of the new fields, one can see that 



the lowest order non-vanishing terms are of the order 
Il^Li r i\a, where the product runs over all channels of a 
given contact: 



Z N ~ Y\ J dr a cos ( ^y")fli,i(Tq)^ /exp ^2m \ 6i,i a (T a ) - i ^Oi,i, a {r N ) + §i,i,a(Ta')\ ^ I 



Taking the averages (. . . }§ with the actions yields E\ N Jl^Li Il tt '>a( T « — T <*') 2 ^ N f° r the product 

of correlators JT Q (. . . )^ . 

"Sffl, |J = = 7r|u; m |0?i , (A3) Thus, the effective action for reads 

J z — ' a 



AT 1 



E, 



H)>f e ft™./*, n 5rz^-~(^EM T .) 



.(A4) 



r 



Note that it is important to keep the non-local in time 
structure of the cosine-term. 

At this stage, we can proceed to integrate out all the 
remaining modes except the static one, 6 m =o ~ as in the 
single-channel case. The prefactor of the cosine-term 
Vq — (E c /ir) Y[ a Ti\a is renormalized according to 

V - V(T) = V Q exp{ - ^ EW T «W r «'))<Wo} 

a.Oi' 

= Vq exp I - ^ f(u m )(l + Jr E cosw mT- QQ ')|, 



a,a'>a 



where = T/(AM 2 ) £ q {(£ q + ^knl)" 1 + 

(Trlwml)" 1 }; see section ITT1 Since typical time differences 
T aa ' = T a — T a r are of the order 1/T (the time integrals 
are dominated by the upper limit of integration) , the last 
cosine-term inside the expon ent m ay be disregarded. As 
a result we find V(T) = V y/T/E c . 

Finally, one may perform the multiple time integra- 
tions in the prefactor of the cosine. The integral over the 
center-of-mass time r — ^ Q r a /N contributes a factor 
1/T, while the integration over TV — 1 independent time 
differences r a — r yields a constant cat multiplied by the 
logarithmic factoid \xiE c /T. The latter follows simply 
from power counting. The same logarithmic factor ap- 
pears in the framework of the phase model, Appendix iBl 
as a result of zero-mode integration. Since all our eval- 
uations of 7 are done up to a numerical factor, we shall 
not keep this logarithm explicitly. We, thus, reproduce 



Eq. 



with 7 (T) ~ ^TjEc I! 



N 



T < To follows the same way as discussed in the main 
text for the single-channel case. 

APPENDIX B: PHASE MODEL 

In this Appendix, we establish correspondence between 
the charge representation, employed in the paper, and 
the more commonly used phase model. The latter may 
be straightforwardly derived starting from the fermionic 
tunneling Hamiltonian. Integration over the fermionic 
degrees of freedom (under the assumption of vanishing 
level spacing in every grain, 8 — > 0), leads to a model for- 
mulated in terms of the dynamic phase variable^, 4>\(t). 
Its time derivative, </>,(t), has the meaning of a fluctuat- 
ing instantaneous voltage on grain 1. The resulting action 
is a straightforward generalization of the Ambegaokar- 
Eckern-Schon (AES) action^ to the array geometry. It 
consists of the charging term, 



SM= Jdr^2 

n 1 



(Bl) 



and the dissipative term 



5 dM = ^/L d r'E Sin2( ^ (r) "^' (T,)) 



<U') 



sin^(7rT(r-r')) 



Continuation to 



(B2) 

describing tunneling between nearest neighbor grains 
(1,1'>. Here, 0ir = (0,-0,0/2. 
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The phase field 0i(r) obeys the boundary condition 
4>i(fl) ~ 01(0) = 2ttWi, where W\ G Z is an integer called 
winding number. In addition to the trivial configuration 
4>\ = 0, the stationary configurations of the dissipative 
action, Sa, are given by Korshunov instantons 24 



nf 



(B3) 



a=l 



characterized by the spatially-dependent winding num- 
ber W\ and a set of complex parameters \z a \ < 1. In 
the regime T -C E c , these configurations are a good ap- 
proximation to the saddle points of the total action. In 
the same approximation the z a are zero-mode coordi- 
nates: the instanton action is almost z-independent (safe 
for the charging terms that weakly depends on z). Ne- 
glecting this dependence, the action for a certain winding 
number configuration, {Wi}, reads 



tt 2 T 



(B4) 



<U'> 



where Ww = W\ — Wy. In the regime, we are inter- 
ested in, T <SC E c and j > 1, the dominant contribution 
comes from the second (tunneling) term in this expres- 
sion. Since the latter depends only on the differences of 
winding numbers on neighboring grains, it favors con- 
figurations with spatially extended regions with a fixed 
constant winding number, e.g. W\ — ±1 for a closed set 
of grains 1. We shall refer to such sets of grains with a 
fixed non-zero winding number as "islands" . A typical 
phase-field configuration contains, therefore, a number 
of "islands" (with fixed non-zero windings) embedded in 
the sea of W = grains. An example of such a config- 
uration is shown in Fig. 0] According to Eq. (|B4(1 , the 
action cost of one such island with the winding number 
W is S = AW 2 {-k 2 T/E c ) + L\W\{g/2), where A is the 
area of the island (number of grains inside) and L is its 
circumference (number of contacts with a winding num- 
ber jump across them). 




FIG. 4: A typical island configuration is shown. The numbers 
correspond to winding numbers in the phase model. 

The same picture was employed in previous studies 
of 2d granular arrays^, where the island structure was 
mapped onto the, so called, solid-on-solid model. The 
latter is known to exhibit a BKT transition. [A true tran- 
sition takes place if there is no cost for the island's area, 



as is the case for the model with mutual capacitances 
only. Indeed, in such a model both charging and tunnel- 
ing terms provide a cost proportional to the island's cir- 
cumference, L. In presence of the self-capacitance (and 
thus an area-proportional cost) the transition is smeared 
into a crossover.] What was missed in the previous stud- 
ies is an account of fluctuations on top of the stationary 
island-like configurations. 

We provide such an account here. Consider a station- 
ary configuration consisting of a single island with a fixed 
winding number W. Expanding to the second order in 
deviations cj)\ — cp^ 1 ' + ipi, one finds for the fluctuating 
part of the action 



m 1,1' 



(W) 
\,V ,m 



(B5) 



where m is a Matsubara index, and My^^L = M^f' is 
the fluctuation matrix. The fluctuation factor associated 
with this configuration is given by 



n 



detM 



(0) 



\ det M 



(W) 



exp -^InfM^r MW 



(B6) 

where is the corresponding fluctuation matrix for 

the flat (W — 0) stationary configuration. In the regime 
g 1, the dominant fluctuation contribution comes from 
the expansion of the tunneling term, S^- This leads 



to M, 



(o) 
i,i', i 



Ah 



(0) 



= — g\u) m \ for nearest neighbors (1,1'), while 

l.i. in - Ei^i^5, ra and M$ >m = otherwise. In 
presence of the island, the off-diagonal elements of the 
fluctuation matrix are changed to M^^ m — — g\u) m _\w\\ 
(and the diagonal elements accordingly), only if 1 and 1' 
are nearest neighbors laying across the island's boundary. 
As a result, one may write Mff' = Mm — | W|<5M, where 
the matrix <5M has entries ±27rTg along the island's 
boundary and zeros everywhere else. Returning to the 
calculation of the fluctuation factor, Eq. I|B6|I . one finds 



Tr In M. 



r (0) 



M 



(W) 



Tr In 



1 



\w\ (m. { ° ] ) W 



— |W|Tr Tm^I SM. Higher order terms in the ex- 
pansion of the logarithm are rapidly convergent upon 
Matsubara summation and, therefore, may be safely ne- 



r(0) 



summation over the Mat- 



glected. Since Mi, 
subara index in J2 m \W\Tr [M^ SM leads to the log- 
arithmic divergence^. It is cut off by the charging part 
of the action at to « gE c /T ^> 1. The summation (trace) 
over spatial indices results in a factor proportional to L, 
the island's circumference, as it counts the number of 
non-zero entries in <5M. Finally, a careful evaluation of 
the numerical coefficient 26 leads to the fluctuation factor, 
Eq. JB6J, equal to (gE c /T) L \ w \/ 2 . 

As a result, an island of winding number W with area 
A and circumference L contributes to the partition func- 
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tion of the model with the relative factor 

(Actually, the statistical weight of an island contains 
also a factor (In E c /Ty w > , coming from the zero-mode, 
z a , integrations^. This factor has its exact analog in 
the charge model, mentioned at the end of Appendix lAl 
Hereafter we omit it for brevity.) 

We shall show now that the pcrturbative expan- 
sion in powers of j(T) of the classical model, Eq. JSJ, 
leads to the same island picture. In this case, ev- 
ery island carries the relative factor Pw(A, L) = 

L-* 2T l E °y W \Ecl(T)/(2-K 2 T)) Lm . We can, thus, 

identify the two models provided j(T) ~ ^JgT/ E c e~ g / 2 . 
Notice that j(T) cx y/T/E c is exactly what one expects 
for the high-temperature, T > To, charge model. At 
lower temperature, non-linear fluctuation corrections in 
the phase model diverge 4 " and the above treatment runs 
out of validity. However, having establish the equivalence 
of the phase and charge models at T > T , one may 
proceed with the analysis of the latter even at smaller 
temperatures. 

To complete the proof, we elucidate now the island 
structure of the perturbative expansion of the charge 
model, Eq. (JSJ). Consider the expansion of the parti- 
tion function Z — J DO exp|— with the action 

S[0] given by Eq. Q, in powers of the small parame- 
ter 7. The partition function can be written as Z — 
EnLo Z n (EcY/T) n , where Z n is a product of n cosine 
terms averaged with the action S c [&\ = E c ^(V • 9\) 2 /T. 
There are two types of contributions: a) terms with 
higher powers of the cosine taken at the same link and b) 
terms involving different links. The first class of terms 
describes pcrturbative corrections to the conductance of 
a single contact, and may be shown to be equivalent to 
those of Ref. [4| in the framework of the phase model. 



The second class corresponds to the instanton terms and 
is the subject of our focus. These terms exhibit "island 
formation" . To illustrate this, let us label the coefficients 
7i,i (even though we assume them all to be equal), where 
i = x,y. Terms of the form F| "/ a are non-zero only 
if the lines crossing all contacts a = (i,l) form closed 
loops, see Fig. 0] I.e. the lowest-order non-local term 
is proportional to 7 4 = J x ,a x ,i +eai lv,l7v,l+e v ~ involving 
all the four links surrounding grain 1. This property of 
the model is due to the presence of massless modes, as is 
explained below. 

Rewriting 9 = Vx + V x 77, one finds that the charg- 
ing action takes the form S c \x, rj\ — E c J2i(^Xi) 2 /T 
and is thus ^-independent. The rotational field rj is 
therefore strictly massless. As a result, as long as 
an argument of the cosine (exponential) function con- 
tains the 77-field, it averages to zero. Indeed, to ob- 
tain Z n one has to average expressions of the form 
exp{2 7 ri(X; j S x= i(±)^,i : , x +E" B = a +i(±)0w,ii v }> where j, 
labels contacts in i-direction. The terms containing r) in 
the argument of this exponent, vanish. Therefore non- 
vanishing are only those terms that have Y) j ^ivij^ ~ 
nh*-e*) + £j B ±{ ^h v ~ ^ijy-ej = 0. It can be seen that 
this condition corresponds to the island structure. As a 
result, every island brings a factor (E^ /T) L , where L 
is its circumference, that simply reflects the order of the 
perturbation theory needed to create the island. For a 
proper (i.e. island-like) term of the perturbation theory, 
the averaging over the massive x~fields results in the fac- 
tor exp{— tt 2 TA/E c }, where A is the area. Finally, the 
integer index \ W\ corresponds to the possibility of having 
a non-zero term of the perturbation theory, where links 
surrounding an island are included \W\ times each. We 
have shown, thus, that the perturbation theory in the 
charge model, Eq. JHJ), produces the same island struc- 
ture as the instanton expansion of the phase model - with 
the same relative factors, Eq. (|B7|I . This completes the 
proof of the equivalence of the two models and provides 
the value of j(T), Eq. @, for g > 1. 
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